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Form factors for local spin operators of the XXZ Heisenberg spin-g finite chain are 
computed. Representation theory of Drinfel'd twists for the quantum affine algebra 
Mq(sh) in finite dimensional modules is used to calculate scalar products of Bethe 
Q\ ■ states (leading to Gaudin formula) and to solve the quantum inverse problem for 

local spin operators in the finite chain. Hence, we obtain the representation of the 
n-spin correlation functions in terms of expectation values (in ferromagnetic reference 
state) of the operator entries of the quantum monodromy matrix satisfying Yang- 
Baxter algebra. This leads to the direct calculation of the form factors of the XXZ 
Heisenberg spin-^ finite chain as determinants of usual functions of the parameters 
of the model. A two-point correlation function for adjacent sites is also derived using 
similar techniques. 
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1 Introduction 



One of the most challenging problems in the theory of low dimensional quantum integrable mod- 
els p], 0, S, 0, |5L Hj, after finding the spectrum and eigenstates of the corresponding Hamiltonians, 
is to construct exact and manageable expressions of their form factors and correlation functions. 
This is a fundamental problem both to enlarge the range of applications of these models in 
the realm of condensed matter physics and to better understand their underlying mathematical 
structures. Until recently, only very few models were known for which correlation functions can 
be computed exactly. Typical examples are the Ising model (related to free fermions) 0, ^, |9| 
and conformal field theories pT[ . 

Beyond these models, in the framework of integrable systems solvable by means of Bethe 
Ansatz [12, 1^, 14, 2], [l], ||, related to a Quantum Group structure [15, 16, [L7], 18] and associated 
to an R-matrix solving the Yang-Baxter equation, one can distinguish at present essentially two 
different but complementary approaches that have been designed to deal with this problem. 



One of them relies on the study of form factors and correlation functions of quantum inte- 
grable models directly in the infinite volume limit. The roots of this approach are twofold: 

On the one hand, it comes from the study of analytic properties and bootstrap equations for 
the factorized S-matrices and form factors of integrable quantum field theories in infinite volume 
fli^ , H|. Typical models here are the two-dimensional Sine-Gordon relativistic quantum field 
theory and the Non-Linear Schrddinger model. There it was realized that the set of equations 
satisfied by the form factors are closely related to the q-deformed Knizhnik-Zamolodchikov 
equations arising from representation theory of quantum affine algebras, and their q-deformed 
vertex operators |(| |2l], ||, ||, |4], |§, gf| |3- 

On the other hand, it uses the Corner Transfer Matrix introduced by Baxter [28, E9j, |], |30| 



in the context of integrable models of statistical mechanics, and the relation of its spectra to 



characters of affine Lie algebras [31]. Typical examples here are the six- vertex model and the 
XXZ Heisenberg spin-^ infinite chain. In such models, using very plausible hypothesis about 
the representation of the Hamiltonian as a central element of the corresponding quantum affine 
algebra (here U q {sl2)) in the infinite volume limit, the space of states is constructed in terms 
of highest weight modules of U q {sl2) [32, 33, 0], the combinatorial aspects of this construction 
being related to the theory of crystal bases [Bj] , 34]. Form factors and correlation functions 
are then described in terms of q-deformed vertex operators, leading via bosonization [35], to 
integral formulas for them. As a result of its algebraic formulation, very parallel to the one in 
conformal field theory, there was a rapid development of this approach, leading to rather explicit 
expressions for correlators and their short distance behaviour Q. 

It should be mentioned however, that the application of this method seems more difficult 
for time or temperature dependent correlators, or, for quantum spin chains, if an external 
magnetic field is present, namely in situations where the infinite symmetry algebra is not as 
clearly identified as in the pure case. 



In these directions, one should cite [36, 
particular a finite volume analysis of form factors has been undertaken. 



where in 



The other approach, described essentially in the book || , is based on the detailed analysis of 
the structure of Bethe eigenstates and in particular of their scalar products properties. One of the 
starting points of this approach is the Algebraic Bethe Ansatz (or Quantum Inverse Scattering) 
method fll4| , |lj || and the derivation in this framework by Korepin of the Gaudin formula for the 
norm of Bethe eigenstates |}q]. Then, to overcome the enormous combinatorial complexity due 
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to the structure of Bethe eigenstates, the two key ingredients of this method are on the one hand 
the so called "dual fields approach" [^] and on the other hand the determinant expression for the 
partition function with "domain wall" boundary conditions [40]. Using these auxiliary quantum 
"dual fields", determinant representations of correlation functions are obtained [41, 42, |43|, H, 
containing however vacuum expectation values of these auxiliary "dual fields" , which cannot be 
eliminated in the final result. Hence, explicit expressions for the correlators cannot be obtained 
directly from this approach. Instead, the strategy is to embed these determinant formulas in 
systems of integrable integro-difference equations from which only large distance asymptotics 
of the correlation functions can be extracted from the resolution of (matrix) Riemann-Hilbert 
problems. 

Let us nevertheless note here, that in simpler models, in particular at so-called free-fermion 
points, such as the XXO model or the Non-Linear Schrodinger model at infinite coupling constant, 
more explicit results can be obtained 44, 45]. 

Although mainly restricted to the determination of these large distance asymptotics of the 
correlation functions, the very general formulation of this method allows one to apply it to a 
large variety of integrable models, and to correlation functions depending on time, temperature 
and eventually, in the case of spin chains, on an external magnetic field. 

A more algebraic understanding of the Bethe Ansatz approach to correlation functions is 
certainly needed to avoid the combinatorial difficulties encountered in this method, in particular 
if one would like to obtain explicit expressions for the correlators, namely, without auxiliary 
"dual fields". 

In this direction one should mention |4(| where Gauss decomposition of operators was used 
within the Gaudin model to produce an explicit determinant formula for the norm of Bethe 
states, or also [26] where the Gaudin formula follows from semi-classical asymptotics of the 
q-deformed Knizhnik-Zamolodchikov equations. 



The present state of the problem is such that, despite the great advances we just briefly 
described, algebraic derivation of form factors and correlation functions in an explicit and man- 
ageable setting, even for the most elementary models such as the XXZ Heisenberg spin-^ finite 
chain, still poses a formidable problem. 

The main motivation of this article is precisely to understand from a more algebraic point of 
view the Bethe Ansatz approach to correlation functions for finite systems, and to try eventually 
to relate the two above approaches in taking the thermodynamic limit. For that purpose, we 
will mainly concentrate on one of the a priori most elementary models in this context, the XXZ 
Heisenberg spin-^ finite chain. We will show how to compute explicit determinant formulas, 
namely in terms of usual functions of the parameters of the model and without any auxiliary 
"dual fields" , for the form factors of local spin operators (i.e. their matrix elements between any 
two Bethe eigenstates) and for the adjacent sites two-point correlation function. In fact, we will 
also obtain this result for the completely inhomogeneous XXZ Heisenberg spin-^ finite chain. 

Our approach to form factors and correlation functions for this model decomposes into three 
main steps: 

i. we compute representations for scalar products of an arbitrary Bethe eigenstate with any 
other state in terms of determinants of elementary functions of the parameters of the 
model. 
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ii. we solve the quantum inverse problem for the completely inhomogeneous XXZ Heisenberg 
spin-^ finite chain, namely, we reconstruct the local spin operators at any site i on the 
chain in terms of the elements of the quantum monodromy matrix of the chain. 

Hi. we combine these two results to obtain determinant formulas for the form factors of the 
local spin operators, and for an adjacent sites two-point correlation function. 

The key ingredient of our method is the article p7| ]. There a factorizing Drinfel'd twist F 
was constructed and studied. This twist is associated to the iV-fold tensor product of spin-i 
(evaluation) representation of the quantum affine algebra U q (sl2) associated to the completely 
inhomogeneous XXZ Heisenberg spin-i chain of length N . It has been shown in particular that 
the change of basis in quantum space of states generated by this twist F is particularly convenient 
to study the structure of Bethe eigenstates in the framework of Algebraic Bethe Ansatz. The 
main explanation of this is certainly the fact that the F-basis determines a completely symmetric 
presentation of the monodromy matrix operator for the (inhomogeneous) chain, such that the 
action of the symmetry group is trivial in this basis. 

As a result, while creation and annihilation operators of Bethe states, namely the operator 
matrix elements of the quantum monodromy matrix, B{\) and C(A), are in the original basis 
represented as huge sums, containing up to 2^ terms, each of them being a product of iV spin 
operators along the chain, their representations in this new F-basis simplify drastically. Indeed, 
in this basis, they are given as sums of only N terms, each of them being simply a local spin 
operator at some site i of the chain, dressed by a pure tensor product of diagonal operators acting 
on the other sites (see section 2 and |47j]). 

This means that the F-basis already solves the combinatorial problem of describing creation 
and annihilation operators of Bethe states. Moreover, we will show in this paper how it also solves 
the combinatorial problem of describing Bethe eigenstates generated by products of creation 
operators -B(A&) on a reference (ferromagnetic) state. This will enable us to compute in an 
explicit way scalar products of a Bethe eigenstate with any other state. The result is obtained for 
the completely inhomogeneous XXZ Heisenberg spin-^ finite chain, as determinants of functions 
of the parameters of the model, and solves the above point (i). 

Point (ii) is also solved by a careful study of the particularly simple structure of the quan- 
tum monodromy matrix in the F-basis. The reconstruction of local spin operators in terms of 
the operator matrix elements of the quantum monodromy matrix is then obtained in a basis- 
independent way. 

Point (iii) uses only the combination of this two results and of some resummation formulas 
we will explain in the main text. 

This article is organized as follows: in section 2 we recall the definition of the inhomogeneous 
XXZ Heisenberg spin-| finite chain, and the algebraic ingredients we will use in the following, 
such as the quantum i?-matrix and the associated quantum monodromy matrix. Further, we 
describe briefly formulas for the factorizing twist F from [47] and its essential properties to be 



used in this article. In particular, we give there the expression of the quantum monodromy 
matrix in the F-basis. In section 3 we derive an explicit formula for the scalar product of an 
arbitrary Bethe state with any other state. Details of the proofs for this section are contained 
in the three appendices at the end of this article. Section 4 is devoted to the solution of the 
quantum inverse problem for the local spin operators. Finally, the main results of our work 
concerning form factors of the local spin operators are presented in section 5. Conclusions and 
perspectives are given in section 6. 
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2 The XXZ Heisenberg spin-| inhomogeneous finite chain 

In this paper, we shall calculate form factors for the Heisenberg XXZ and XXX spin-^ chains 
of length N. The XXZ Heisenberg model is given by the following Hamiltonian: 

N 

HXXZ = J ^2 fm a ™+l + tt+l + A ( Cr m <T m+l ~ !)}> I 2 - 1 ) 
m=l 

the particular case A = 1 corresponding to the XXX chain. We impose here periodic boundary 
conditions. 

Our method is based on the Algebraic Bethe Ansatz [fTi], ^], the central object of which is 
the quantum i?-matrix. For the XXX and XXZ models it is of the form 



R(\,H) 



(I 0\ 

b{\,fi) c(A,/x) 

c(A,m) 

\0 1/ 



(2.2) 



where 



b(X,,)= ^-'K , (2-3) 



<p(\- li + rj)'- 
<p(\-H + r))' 



c(A,m)= ,/ (7?) , „ (2-4) 



with the function (/? defined as 



(p(X) = A in the XXX case, (2.5) 

<p(\) = sinh(A) in the XXZ case. (2-6) 

The i?-matrix is a linear operator in the tensor product of two two-dimensional linear spaces 
V\ <8 V2, where each V{ is isomorphic to C 2 , and depends generically on two spectral parameters 
Ai and A2 associated to these two vector spaces. It is denoted by _Ri 2 (Ai, A2). Such an i?-matrix 
satisfies the Yang-Baxter equation, 

-Rl2(Ai,A2) -Rl3(Al,A3) i?23(A2,A3) =i?23(A2,A3) i?13(Al,A3) #l 2 (Ai,A 2 ), (2.7) 

the unitary condition (provided 6(Ai, A 2 ) 7^ ±c(Ai, A2)), 

i?i 2 (Ai,A 2 ) J R 21 (A 2 ,A 1 ) = 1, (2.8) 

and the crossing symmetry relation, 

(7 ® 1) i? 12 (Af, A 2 ) (7 ® 1) = 4 1 i(A 2 , Ai) p(Ai, A 2 ), (2.9) 

with /9(Ai,A 2 ) being a scalar function, and 7 a 2 x 2 matrix such that 7 2 = 1, 7* = ±7, the 
upperscript t,- meaning the usual transposition of matrices in the corresponding space (j). For 
the rational case, 

Af = A x — 77, 7 = cr y , p(Ai,A 2 ) = A \ A2 - , 

Ai — A 2 
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and for the trigonometric case, 

sinh(Ai — A2 — rj) 
sinh(Ai — A2) 



A? = Ai - r? + i7r, 7 = cr x , p(Ai,A 2 ) 



where a x and a y are the standard Pauli matrices. 

Identifying one of the two linear spaces in the i?-matrix with the two-dimensional Hilbert 
space Tin of SU(2) spin-^ corresponding to the site n of the chain, it is possible to construct 
the quantum L-operator of the model at site n as 

L n (A,£n) = R n(X,( n ), (2.10) 

where £ n is an arbitrary (inhomogeneity) parameter dependent on the site n. The subscripts 
mean here that R$ n acts on the tensor product C 2 <g> H n . The quantum monodromy matrix of 
the total chain defined as the ordered product of L-operators is given by 

MX) = T 0tUM (\; t N )=R 0N (\,Z N )...R ox (\,Z 1 ). (2.11) 

It can be represented in the first space as a 2 x 2 matrix, 



(A{X) B{X) 
{ ' \C{X) D{X) 



(2.12) 



whose matrix elements A(X) = Ai...n(\;£i, . . . , £jv), B(X) = 5i..jv(A; £i, . . . , £at), C(X) = 

Ci...n(X; Ci, ■ ■ ■ , £,n), D(X) = Di...n(X; . . . , £jv) are linear operators on the quantum space of 

N 

states of the chain Ti = ® TC n . Their commutation relations are given by the following relation 

n=l 

on C 2 ® C 2 : 

R 12 (X,fi) Ti(A) T 2 (ji) = r a ( M ) Ti(A) i2 12 (A,//), (2.13) 

with the usual tensor notations T\(X) = T(A) <8> Id and T 2 (/i) = Id <S> T(fi). 

The monodromy matrix satisfies moreover the following (crossing symmetry) relation leading 
to the definition of the quantum determinant: 

T ,i..jv(A;£i,... ,6v) 7 o T*° ^ N (X S ; 6, . . . ,£n) 7o = p(A;6,--- ,£n) 1, (2.14) 

where p(A;£i, . . . ,£n) = YiiLi p{X,£i), and 70 is the matrix 7 of eq. ( |2.9D acting in space 0. 

One also defines the transfer matrix T(A) as the trace A(X) + D(X) of the total monodromy 
matrix. Thanks to the Yang-Baxter equation and the invertibility of the i?-matrix, the transfer 
matrices commute with each other for different values of the spectral parameter A. For the 
homogeneous case where all parameters £j are equal, the Hamiltonian ([O]) can be obtained in 
terms of the transfer matrix by means of trace identities. 

The Algebraic Bethe Ansatz, which deals with the problem of diagonalizing simultaneously 
T(A) for all values of A, supposes the existence of a reference state | 0), called pseudo- vacuum, 
such that 

A(A)|0)=o(A)|0), 
£>(A)|0)=d(A)|0>, 
C(A)|0) = 0, 
B(X)\0}^0. 
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For the XXX or XXZ model, the pseudo-vacuum is the completely ferromagnetic state with all 
the spins up, and a(X) = 1, d(X) = YliLi Common eigenstates of the transfer matrices 

for different values of the spectral parameter A are obtained as successive actions of operators 
B on the pseudo-vacuum Y!j=i 0); f° r an Y se ^ of n spectral parameters {Xj, 1 < j < n} 

solution of Bethe equations 

n = x ' i -*- n> (2 - i6) 



3 

3^k 



with 



The corresponding eigenvalue for the transfer matrix is then 

n n 

t(jm, {Xj}) = a{n) [] b~ 1 (Xj,fJ,) + d(jj,) ]J b' 1 ^, Xj). (2.18) 

3=1 3=1 



Let us now turn to the description of the key object we will use to compute scalar products 
of Bethe states and to solve the Quantum Inverse Problem for local spins, leading finally to the 
form factors formulas: the factorizing F-matrix associated to the above i?-matrix. 

The concept of factorizing F-matrices was defined in E7|, following the concept of twists 



introduced by Drinfel'd in the theory of Quantum Groups [15]. To be essentially self-contained 
we briefly recall here their main properties and refer to ]^7| for more details and proofs. 

Due to the Yang-Baxter equation and to the unitarity of the i?-matrix associated to the 
XXX and XXZ models, for any integer n one can associate to any element a of the symmetric 
group S n of order n, a unique i?-matrix -Ri... n (£i, • • • , £n) constructed as some ordered product 
(depending on a) of the elementary i?-matrices Rij{£,i,£,j) (see ^7|), such that 

^l...n(£l> • • • i£n) T 0>1 _ n (X; £l, ... ,Cn) = 7o,<T(l)...<r(n) £<r(l) > • ■ ■ ,£a(n)) R l...n(Cl, ■ ■ ■ ,Cn)- 

(2.19) 

A factorizing F-matrix associated to a given elementary R matrix is an invertible matrix 
Fi...n(£,i, ■ ■ ■ >£n) satisfying the following relation for any element a of S n : 

Fa(l)...a(n)(£,a(l), ■ ■ ■ , £<r(n.) ) ^l...n(£l> • • ■ , £n) = F l-n (6 , • • • ,£n)- (2.20) 

In other words, such an i^-matrix factories the corresponding i?-matrix. Taking into account 
the fact that the parameters £j are in one to one correspondence with the vector spaces Vi, we 
can adopt simplified notations such that 

-Fl...n(£lj • • • > £n) = Fl...m 
Fo-(l)...o-(n)(t,a(l), ■ ■ ■ ,£a(ri)) = F a(l)...a{n) j 
Fl,2...n(£,l',£,2 ■ ■ ■ , £n) = Fl,2...n- 
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An explicit formula for a triangular F-matrix corresponding to the XXZ model has been con- 
structed in [JiTfl. It reads for any integer n, 

Fi... n = F 2 ... n F h2 ... n , (2.21) 

= Fn-ln F n _2,n-ln ■ ■ ■ ^1,23.. .ni (2.22) 

where the partial F-matrices i ? i,j+i...n(Cij £i+ij • • • >£n) are given in terms of the /^-matrices as 

(11) (22) 

Here we have defined the partial i?-matrices acting vnVi® ■ ■ ■ ®V n as 

■fti,i+l...n(£i)£i+ls • • • ,£n) = Rin(£ii Cn) ■ • • (2.24) 

and is the elementary matrix e^ fc ^ acting in space i, with matrix elements e^j = 5 a k&u- 
The partial F-matrix i*b,i... tv(A; £1, • • • , £jv) has a useful expression as a 2 x 2 matrix in the first 
space in terms of elements of the quantum monodromy matrix: 

iV...Ar(A;£i, • • • ,£n) = ( r c \ r> f\-e c \ I ' (2-25) 

V ^l...AT(A,4i, . . . ,4at) Dl...N\A,Sl> ■ ■ ■ >4iVj /r i 

Let us note here two important properties we will use in the following. The first one can be 
derived directly from the above relations between R, T and F, leading to 

Fi...n(£i, ■ ■ ■ ,£n) r 0) i...jv(A;^i, . . . ,£jv) Fi.„ N (£i, ■ • • >€n) = 

= ^7(l)...cr(Af)(C CT (l) 5 • • • ,£a(N)) T 0,a(l)...a(N)(^,£a(l), ■ ■ ■ ,tcr(N)) F a{l)...a{N) (^(1), • • • ,£o-(JV))- 

Hence, it means that in the i^-basis, the monodromy matrix T defined as 

?o,i...jv(A;£i, . . . ,£jv) = i*i..jv(£i., ■ • ■ ; £iv)2o,i..jv(A; £1, • • • i£n) F^ N (£\, . . . >£zv), 

(2.26) 

is totally symmetric under any simultaneous permutations of the lattice sites i and of the cor- 
responding inhomogeneity parameters 

The second property, proved in [|47[ , is as follows: for the XXZ-i model, the quantum 
monodromy operator is a 2 x 2 matrix with entries A, B, C, D which are obtained as sums of 
2 _1 operators which themselves are products of N local operators on the quantum chain. As 
an example, the B operator is given as 

N 

#i..jv(A) = crT" ttj + crT" (aj a^) ftijk + higher terms, (2.27) 
i=l i^j^fe 

where <r + , cr~ and sgz are the standard Pauli matrices and the matrices f2j, are diagonal 
operators acting respectively on all sites but i, on all sites but i, j, k, and the higher order terms 
involve more and more exchange spin terms like a~ a^. It means that the B operator returns 
one spin somewhere on the chain, this operation being however dressed non-locally and with 
non-diagonal operators by multiple exchange terms of the type a~ at ■ 
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So, whereas these formulas in the original basis are quite involved and cannot be used in 
direct computations, their expressions in the F-basis simplify drastically. From |^7j we have 

Al..jv(A;£i, • • • ,6v) = Fi_ N (£i, . . . ,£ N ) F>i..jv(A;£i, • • • ,€n) -FfV^i 
The operator B representing the operator B in the F-basis is given by 

Similarly we have for the operator C, 

i=i ^ ^ u 1 / b1 

and the operator A can be obtained from quantum determinant relations ( |2.14| ). 

We wish first to stress that the operators A, B, C, D satisfy the same quadratic commutation 
relations as A, B, C, D. Second, each of the operators B and C is reduced to an elementary 
sum on the sites of the chain of the corresponding spin operator at each site dressed diagonally, 
which is to be compared to their expressions in the original basis where they are given as sums 
of 2 N terms involving much more complicated operators. 

It really means that the factorizing F-matrices we have constructed solve the combinatorial 
problem induced by the non-trivial action of the permutation group Sn given by the i?-matrix. 
In the -F-basis the action of the permutation group on the operators A, B, C, D is trivial. 
Moreover the operator A + D which contains the Hamiltonian of the model together with the 
series of conserved quantities, is now a quasi-bi-local operator. 

Further, it can be shown that the pseudo-vacuum state is left invariant, namely, it is an 
eigenvector of the total F-matrix with eigenvalue 1. Hence, in particular, the Algebraic Bethe 
Ansatz can be carried out also in the F-basis. For the scalar products of the quantum states of 
the model we have 

( | C(Ai) . . . C(X n ) B(X n+1 ) . . . B(\ 2n ) | ) = 

= <0| C(X 1 )...C(X n ) B(\ n+l )...B(\ 2n ) |0). (2.31) 

Hence, thanks to these very simple expressions, a direct computation of Bethe eigenstates and 
of their scalar products in this F-basis is made possible, while it was completely hopeless in the 
original basis. There, only commutation relations between the operators A, B, C, D can be 
used, leading (see ||) to very intricate sums over partitions. 

We now end this section with some useful formulas making the computation of the F-matrices 
simpler, namely the expressions of the partial F-matrices in the F-basis. This will help us to 
solve the quantum inverse problem for the local spin operators. 

Factorizing F-matrices are given in (2.22) as an ordered product of partial F-matrices like 
Fi 5 2...n- These object are constructed in terms of the i?-matrix i?i 5 2...n- However this quantity is 



> ••• (Cat) 

(2.28) 
(2.29) 
(2.30) 
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highly non trivial to compute explicitly, since it involves in fact sums of 2 n_1 terms. In contrast, 
the partial F-matrices in the -F-basis can be obtained explicitly, while they also lead to the 
construction of factorizing F-matrices F\2...n- We have (using again simplified notations) 

Fl„.n = Fl,2...n F2...T1 

= F\,2...n F2,3...n • • • -Pra-ln, (2.32) 

and the partial F-matrix F\2...n reads as a 2 x 2 matrix in the first space 1: 

^l,2...n(£i;£2, • • • ,£n) = F 2 ... n (£,2, ■■■ ,€n) ^l,2...n(6; 6; • • • > £n) -^SUnfe, ' " > £") 



1 

C2...n(Ci;^2, • • • ,£n) £>2...n (6 ! 6 , • • • , £n 



(2.33) 
(2.34) 



[i] 



It is a very simple object to compute from the formulas of this section. Hence using the F-basis 
we have also obtained a more explicit and elementary formula for the F-matrix itself. 

3 Scalar products of Bethe states and the Gaudin formula 

In this section we calculate the following scalar products of states constructed by the action of 
the operators B(X) on the pseudo- vacuum, 

n n 

S n ({ f i J },{X k }) = (0\ IJCCMi) 1°)' C 3 ' 1 ) 

j=l k=l 

when one of the sets of parameters, for example {A^}, is a solution of Bethe equations. Hence 
the state n&=i B(X k )\ 0) is supposed to be an eigenvector of the transfer matrix, 

n n 

{A{n) + J] B ( X k)\ ) = t(h, {\ k }) B(\ k )\ 0), (3.2) 

fc=i fc=i 

with the eigenvalue 

n n 

r(fi,{X k }) = a(^Hb- 1 (X k ,fi) + d(fi)l{b- 1 (^X k ). (3.3) 

k=l j=l 

We will prove the following theorem: 

Theorem 3.1. Let {X\, . . . , A n } be a solution of Bethe equations 

and {/xi, . . . ,fj, n } be an arbitrary set of parameters. Then the scalar product ( |3.1| ) can be repre- 
sented as a ratio of two determinants 
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of the following n x n matrices T and V: 





Fab = ^- T 0b, {Afc}), V ab = — — , 1 < a, b < n. (3.5) 



Proof — Let us first note that the computation of the derivatives in the matrix T and of the 



determinant of the matrix V gives a formula obtained in 48]. However, the proofs proposed 



in [48, 49 1 are quite complicated and use some recursion relations for the scalar product or the 
dual field representation. Here we give a direct proof of this formula for XXX and XXZ models. 

The usual approach to the scalar product developed in |38], £il|, Q is based on the commu- 
tation relations between the matrix elements of the monodromy matrix (operators -A(A), B(X), 
C(X) and D(X)). It leads to the recursion relations for the scalar products. Instead of it we use 
the explicit representations for these operators in the i^-basis. Indeed, as the vacuum vector is 
invariant under the action of the operator F, the scalar product ( |3.1| ) can be rewritten in terms 
of the operators in the -F-basis, 

n n 

S n = (0\ \{C{ H ) \{B(X k ) |0). (3.6) 

3=1 k=l 

To perform the computation, it is convenient first to change the normalization of the oper- 
ators -B(A) and C(A): 

We thus want to calculate the "renormalized" scalar product in the F-basis 

S n = (0| C(// n )...C( A t 1 )B(A 1 )...B(A n ) |0), (3.8) 

in which we suppose {A&} to be a solution of Bethe equations. 

The idea is to insert in the scalar product complete sets of states | i\ , . . . , i m ) beyond each 
operator C (A) , where we denote by \i\, . . . ,i m ) the state with m spins down in the sites 
ii, . . . , i m and with N—m spins up in the other sites. We are thus led to consider the intermediate 
functions 

G (m) ({A fc },/ii, . . . ,n m ,i m+ i, ... ,i n ) = (Wli ■■-,in \ C (flm) . . . C (m) 1 (Ai) . . .B (A n ) I 0), 

(3.9) 

the last one being the scalar product, 

G {n) ({A fc },/"l,-..,Ain) =§n- 

There is actually a very simple recursion relation between these functions, 
G (m) ({A fc } ) X 

x G (m ~ 1) ({A fc },/ii, . . . ,/i m _i,j,i m+ i, . . .,i n ), (3.10) 
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where the matrix elements of the operator C (/u m ) can be easily calculated in the .F-basis: 
|C( M )| 

^) l=m+l (3.11) 

The function G , defined as 

n 

G(°\{\ k },h,...,i n ) = (h,...,i n \ Y[»(\ k ) |0>, 

fc=i 

is closely related to the partition function of the six-vertex model with domain wall boundary 
conditions, which was initially given in p0[ . In Appendix A we compute directly this partition 
function using the -F-basis representation for the operators B(X) and C(A). The function G (°) 
is then calculated in Appendix B: 

n n 

n n viK - ii h + v) 

G (°) ({A Q }, . . . , i n ) = a 7l k -\ v n 7T jpr det AA({A Q }, }) , 

where the n x n matrix Af({\ a }, {£i-}) is defined by 

Ma3 = -7i 1 a TT' l<a.i<«- (3-13) 

One should now perform the summation in the relation (3.10) and compute successively the 
functions G( m ). It can be done using some identities for rational functions. Detailed calculations 
are given in Appendix C. One finally obtains, both in the XXX and XXZ cases, 

s a \ s\ w (r{n)<s\ \ \ detiJ({A Q },{/i i }) 

»n({^|,{A Q }) = G { >({X a },fli, . . . ,H n ) - 



n f(^ k - ^) n <p(*0 - a q ) ' ( o u) 

j>k a<0 \ ' ' 

where the matrix elements of the n x n matrix H({X a }, {^j}) are 

Hab = (T~T — ) ( r ^ fo ) II ^( Am ~~ + V) ~ II ^(Am - Mb - 7 7)) • (3-15) 



This formula was originally obtained in 
of eigenvalues of the transfer matrix (|3T3 ) 



48| . It can be rewritten in a very simple form in terms 



n n 



n n <p{pj - A «) 

S n (M,{*c,}) = rr 7~~~~~\ n a --detr({^},{A a }), (3.16) 

n ^(^fc - Mi) n v(A/s - A a ) 

i>fe a</3 

where the matrix T is a Jacobian, 

d 

Tab = -g^T(fl b ,{\ a }). 
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Using a well known formula, 



n <p(k - Xb) n ^(vk - 

detV= a<b - - > <k , (3.17) 

n n ^(m* - x a) 

k=la=l 

with the matrix V defined by ([T^), one can express the coefficient in ( |3.16| ) as a determinant 
and obtain finally the representation ( |3.4[) . 

One should also mention that one can suppose from the beginning ( | []" =1 C(A a ) to be 
a Bethe state (instead of n«=i B(\ a )\ 0)), make almost the same calculations, and obtain the 
following result: 

Sn({A Q }, M) = {Aa})- (3.18) 

□ 



It can be easily seen that taking the limit /j, a — ► A , a = 1, ... , n, in the expression ( 3.14| ) 



one obtains a very nice proof of the Gaudin formula for the square of the norm of the Bethe 
wave function, initially proved by Korepin pSj], 

n n 

n„=(oi n c (^-) n^( x k) io> , 

j=i fe=i 

11 y(A a — A/3) 
where is a n x n matrix the elements of which are given by 



4 Quantum inverse scattering problem for local spins and cor- 
relation functions 

In the previous section we have calculated the scalar products of Bethe states using the F-basis. 
Our purpose is now to compute form factors of local spin operators. One possibility is to write 
them in the -F-basis, which allows us to perform their calculation in the same way as for the 
scalar products. There is however a more simple and instructive way to proceed, which consists 
in solving the quantum inverse scattering problem for the local spin operators, that is expressing 
them only in terms of the operator entries A, B, C, D of the quantum monodromy matrix of 
the model. 



4.1 Local spin operators in the F-basis and quantum inverse problem 

As we just said, the first possible way to compute the form factors is to express the local spin 
operators in the F-basis. Note that it can be directly done for spin operators at the first or the 
last site of the chain. For example, due to the recursion relation Fi...jy = F\^...n F2...N and to 
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the very simple form (2.34) of Fx 2..jVj a direct computation of a product of 2 x 2 matrices in 
the space 1 gives the value for in the F-basis, 

Fi... N (Ti F- X N = 5i„jv(fi;£i, ...,£ N ) oy. (4.1) 

Similar expressions hold for ajj, a\ and o~ z N in the F-basis. 

Thus, taking into account that D is totally symmetric, for a given site i of the chain this 
result can be simply translated into the following formula: 

v:...- 1 ^ CV.i-l = Di-Mti) ( 4 - 2 ) 

where we used the short notation I?i..jv(A) = D\.. jv(A; £1, . . . , £jv) for A = £j. Hence, to 
calculate the operator in the F-basis Fl.jv Cj 7 F^ N , one should evaluate the product of two 
"permuted" F-matrices Fl.jv F~ l Nl It can be considered as the expression in the F-basis 
of some propagator Fr N1 i _ 1 Fi...at, for which the following result holds: 

Lemma 1. Let U\ be the propagator F^ 1 N1 - x Fl.jv from site 1 to site i of the chain. It can 
be written into the two following forms: 

U\ = Ri-l,i...Nl...i-2 ■ ■ ■ R2,3...N1 R\,2...N, (4.3) 

i-1 

= H(Ai.MU + Di...N^a))- (4-4) 

a=l 

Proof — Equality (4.3) comes from the factorizing property for the F-matrix ( |2.20| ) for a cyclic 
permutation a [47]: 

F a ...Nl...a-l Ra-l,a...Nl...a-2 = i*a-l..JVl...a-2> (4.5) 

Equality (4.4) follows from the identity 

A\...n{(,o) + Dl...N{t,a) = Ra,a+l...Nl...a-l(£,a',£,a+l, ■ ■ ■ ,£a-l)- (4-6) 



The proof of ( |4.6|) is based on the remark that Ro a (£ a ) is nothing but the permutation matrix 
Po a of the spaces and a. Writing Ai..jv(£a) + ■Di..jv(£ Q! ) as a trace in the auxiliary space Vb, 
and making Po a act on every factor, we obtain, thanks to the cyclicity of the trace, 

^4-1.. jv(£a) + ^ l..J\r(£a) = tr (R()N(£,a) ■ ■ ■ Roa+l(£,a) P()a R()a-l(£,a) ■ ■ ■ -R(Vl(£a)) > 

= Raa-l\£a) ■ ■ ■ Ral(Ca) RaN{Ca) ■ ■ ■ Raa+l(£a), 
= Ra,a+l...Nl...a-l(Ca : , £a+lj • • • i £a-l)j 

which ends the proof of lemma |l[ Note that in (4.3), (4.4), all factors commute with each others. 
□ 

Remark 4.1. The propagator Ul through the whole chain being the identity, let us notice that 

N 

J] (Al-MZa) + Dl...N(£*)) = R N,l...N ■ ■ ■ Rl,2...N = 1- (4.7) 
a=l 

Thus, the inverse (U{) 1 of the propagator on a part of the chain is nothing but the propagator 
Ui+i = n^=«+i (A + D) (£ Q ) on the remaining part. 
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Remark 4.2. The action of the propagator consists in shifting the beginning of the chain from 
site 1 to site i. For an operator entry X\__.n of the monodromy matrix (X = A, B , C or D), it 
means that 

U{ Xl..jv = Xj...jvi...i-i U{, (4.8) 

which, in terms of the monodromy matrix, can be written 

U{ T 0j i...at = 7o i i...jvi...i_i U{. (4.9) 

The lemma [l] allows us to obtain the value of any local spin operator in the F-basis. For 
example, o~~[ for a given site i of the chain becomes 

i-1 i-1 

F U . N a~ F-' N = H(A + D)(Z a ) • a r ■ J] (2 + D)-\i a ). (4.10) 

a=l a=l 

Using this expression — and similar ones for af , of — it is possible to compute directly the 
corresponding form factor by the same method as for scalar products. 

Moreover, as the expressions of the operators B{X) and C(A) in the -F-basis are quasilocal in 
terms of the o~[ or o~f, it is a way to solve the quantum inverse problem for local spin operators. 
Indeed, a direct calculus in the -F-basis gives the identity 

5(6) ar (A + D)(£ i ) = B&), (4.11) 

which, with ( 4. ICQ , leads to a new reconstruction of . There exists a straightforward proof of 
this last formula, that we expose in the next paragraph. 

4.2 Reconstruction of local spin operators 

In this paragraph, we prove an important result concerning the reconstruction of any local spin 
operator in the inhomogeneous spin chain in terms of elements of the monodromy matrix. 

Theorem 4.1. Local spin operators at a given site i of the inhomogeneous XXX or XXZ Heisen- 
berg chain are given by 

i-1 N 

c-- = H(A + D)(Z a ) ■ B&) • J] (A + D)(£ a ), (4.12) 

a=l a=i+l 
i-1 N 

<rt = Il(A + D)(Z a ) ■ C(&) • J] (A + D)(£ a ), (4.13) 

a=l o=i+l 

i-1 N 

af = H(A + D)^ a ) ■ (A -D)^) ■ J] (A + D)(£ a ). (4.14) 

a=l a=i+l 

The proof of this theorem is a straightforward consequence of the following lemma when xq 
is respectively equal to , o-q and ctq : 

Lemma 2. Let Xi be an operator acting on the quantum spin space Vi. We note xq the corre- 
sponding 2x2 matrix acting on the auxiliary space V$. They are related by the identity 

i — 1 i 

tr (x R 0A ... N (Zi)) = l[(A + D)- 1 (Z a ) ■ x t ■ J] (A + £>)(&), (4.15) 

a=l a=l 

where the trace in the left hand side is taken on the matrix acting in Vq. 
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Proof — Arguments used to prove the case i = 1 are quite similar to those of lemma [T|: 

tr (x #o,i...iv(£i)) = tr (x i?oiv(Ci) • • • #02(6) foi), 
= xi i?uv(Cl) • • -^12(6), 

= Xi Rl,2...N((l), 

and we conclude with lemma [l|. 

To prove the general case, let us notice that in the F-basis i2o,i..jv(£i) is completely symmetric 
into the spaces 1, . . . ,N, which enables us to consider that the chain begins with site i: 

tr (x Ro,i...n(&)) = ^iT.jv tr o(£o Ro,i. ..N^i)) Fi... N , 

= F 1...N tr o( x O Ro,i...Nl...i-l(€i)) Fl...N, 

= F± N Fi-NX...i-l tro(xo Ro,i...Nl...i-l(£,i)) ^i7.7Vl...i-l F1...N, 

= f i7.jv Fi...m...i-i Xi {A + D) (&) /•';...'* , Fi...jv. 

The value of the propagator F~ 1 N1 i _ 1 i*i...jv is given by lemma [I], which concludes the proof 
of lemma 0. □ 



4.3 General formula for correlation functions 

These results make it possible to write a general formula for any /c-point spin-spin correla- 
tion function between two Bethe states for the inhomogeneous XXX- \ or XXZ-i Heisenberg 
chain. Indeed, for any integer k and any subset {ij}i<j<k of {1, . . . , N}, with the convention 
i\ < %i . . . < ik, the correlation function for spins at sites ii, ■ ■ ■ ,ik between two Bethe states 
( I C(fii) . . . C(p ni ) and -B(Ai) . . . B{\ n2 ) | } can be written into the following form: 

( I C( M1 ) . . . C0u m ) all < 2 • • • < B ( X ^ • • • B ^ 2 ) I ) = 

i\ — 1 ri\ N 112 

a=l j=l a =i k +lj=l 

12-1 

x <0| C( Pl )...C( Pni ) ■ • n ( A + D )(^) ■ X%& 2 ) ... 

a=i\+l 

ife — 1 

... J] {A + D)(i a ) ■ X%($ ik ) ■ B(X 1 )...B(\ n2 ) |0>, (4.16) 

a=ife_i+l 

where €j, 1 < j < k, takes the values +, — , or z, X e i being equal respectively to C, B and 
A-D. 



Hence we have reduced the problem of computing any correlation function of the XXZ model 
to a simpler problem written only in terms of the operator entries of the quantum monodromy 
matrix of the chain. 



In the next section, we shall compute explicitly the form factors (k = 1) and the two-point 
correlation functions at adjacent sites (k = 2 and i\ = ii — 1). 
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5 Form factors 



We derive here explicit expressions for the form factors of the local spin operators for the finite 
inhomogeneous XXX and XXZ chains. More precisely, we calculate the matrix elements of the 
operators cr+, cr~ and a z m between two Bethe eigenstates. We also give an expression for the 
simplest correlation function of two spin operators at adjacent sites. 



5.1 Operators a m and cr+ 

We begin with the calculation of the following one-point functions, 

n+l n 

F-(m,{ H },{X k }) = {0\ l[C( H )a- JJ B(A fc ) |0), (5.1) 

j=i k=i 

and 

n n+l 

F+(m,{X k },{^}) = (0\ l[C(\ k )a+ X{B{^) |0>, (5.2) 

k=l j=l 

where {Xk}n an d are solutions of Bethe equations. Using the results of the previous 

sections we prove here that they admit the following representations: 

Proposition 5.1. For two Bethe states with spectral parameters {\k}n an d {/Ujln+i? the matrix 
element of the operator cr~ can be represented as a determinant, 



F n (m, {fij}, {A fc }) 



(pm-ldflj}) 3=1 



n+l 



k=l 

x n n ^-K) d ^ lH ' {m - { "' }Ah}l (5 - 3) 

n+l>fc>j>l l<(3<a<n 



where the coefficients 4> m {{^k}) are 



.({A fc })=nn 6_i ( A *'^)' ^ 
k=ij=i 



and the (n + l) x (n + l) matrix H is defined as 

/ \ / n+l n+l \ 

H ab( m ) = 7 x n ( a ( X b) TT <p(m -h + v)~ d(X b ) TT vifij -h-v)) for b< n + l, 

ma-h)\ fji J (55) 

H ^ {m) = ^ a - U +t^a ~ U ' (5 ' 6) 

The matrix element F+(m, {Xk}, {l^j}) of the operator (7+ admits a similar representation, 

F+(m,{A fe },{^}) = };{Afc}) . (5 . 7) 
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Note that in the homogeneous limit £j = 0, j = 1,...,N, the coefficients 4> m ({\k}) are 
expressed in terms of the total momentum P of the state parameterized by {A^}, 

4>m({^k}) = exp{-iPm}, 

with 

k=l 

Proof - The proof of these representations is rather straightforward. As it was shown in 
previous section (theorem the local operator <r~ can be expressed in terms of the transfer 
matrix and the operator B(£ m ) as 

m-l N 

j=l j=m+l 

Since the Bethe states are eigenstates of the transfer matrix, 

MM + Dm n B(\ k )\o) = ( n b-\\a, & ) n °>> ^ 

k=l \=1 ' k=l 

the product of the operators A(£j) + D(£j) contributes to the function F~(m, {fj.j}, {A&}) as a 
global factor: 

n+1 n 

F-(m,{^},{X k }) = ^ m \{Xk})^ m - 1 ({ H }) <0| YlC( H ) B(U)UB(X k ) |0). 

j=l k=l (5.9) 

Here we used a simple property of the solutions of Bethe equations, 

n N 

nii 6 ~ i ( A *'^)=i- 

k=ij=i 

The right hand side of ( |5.9[ ) thus reduces to a scalar product, 

F~(m, {//,•}, {A fe }) = 4>w{{X k }) <t> m -i{{^ 3 }) S n+ i({fij}, {£ m , Ai, . . . , A n }), 

(5.10) 



which, {/ij} being a solution of Bethe equations, can be computed by means of theorem pA\ 
Hence 

ei-a™ r i r\ i\ i-i/r\ i \ i /r i\ det T n+ i({/Zj}, {£ m , Ai, . . . , A n }) 
(™;{M>l A fc}) = <P m ({Afc})0 m -l({Mj}) 



detV n+ i({nj}, {Cm, Ai, . . . , A n }) ' ^ ^ 



where T and V are (n + 1) x (n + 1) matrices defined similarly as in theorem 3.1 . Writing them 
explicitly one obtains the representation (5.3). 

The form factor F+(m, {A^}, {f-ij}) can be calculated analogously using the representation 
for the operator a m given by theorem [4.1[ □ 
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5.2 Operator a z m 

We calculate here the matrix elements of the operator between two Bethe states, 

n n 

FZ(m, fa}, {X k }) = (0 1 l[C( H )< n^(^)|0}. 

j=i k=i 

We prove the following representation for this one-point function: 

Proposition 5.2. For two Bethe states with sets of spectral parameters {X k } n and {fJ-j} n , the 
matrix element of the operator can be represented as a determinant, 

Ai( m > Wh {Xk}) - 7TVT\ 11 ~T\ c i — V x 

0m-l({Afe} *± tp{\j -im + rf) 



1 



n pink - tij) n <p&p - x c 

j>k a</3 



det n (H({ H },{\ k })-2P(m,{ f i j },{\ k })), (5.12) 



where H({/ij}, {X k }) is the same matrix as for the determinant representation fl3.14|) of the 
scalar product, 



<Mjf) 
(J>a ~ 

and P(m, {Afc}) is a matrix of rank one, 



Hab = ( ^ Z. \ ) ( a ( A &) II ~h + v)- d(Xb) <p(ltj - X b ~ v)j 1 <a,b<n, 



n ^(Afe - a 6 + 77) 

P ab (m) = <p( V ) k=1 — — 1 < a, b < n. (5.13) 

Proof — It is a bit more complicated than for cr _ or o~ + . We use here a representation for 



which follows directly from remark 4.1 and theorem 4.1 



m-l N 
a* m = 2 H(A + D)^) ■ A(U) ■ II ( A + D )^)~^ ( 5 - 14 ) 

j=l j=m+l 

which leads to 

F*(m, {X k }) = 2 ^({Afc}) <f>m-i(M) P ^™, {/%}, {X k }) - S({^}, {X k }), 

with 

n n 

Pi(U,{Hh{Xk}) = (0\ \{C{ H )A{U) l[B(X k )\0). (5.15) 

j=i k=i 

Let us now compute the function P\. It can be done using a well-known formula for the action 
of the operator A(£) on an arbitrary state: 

n n n 

A(U)l[B(X k )\0) = Ylb-\X k ^ m ) JjB(A fc )|0) - 

k=l k=l k=l 
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Hence Pi reduces to a sum of scalar products, therefore to a sum of determinants according to 



theorem It can be rewritten as a single determinant by means of the following formula for 
the determinant of the sum of two matrices one of which being of rank one. Indeed, if A is an 
arbitrary n x n matrix and B a rank one n x n matrix, the determinant of the sum A + B is: 

n 

det{A + B) = det A + det A {j) , 

3=1 

where 

A$ = A ab for b + j, 
A {j) - B ■ 

Using this formula and also the orthogonality of two different Bethe states, one obtains the 
determinant representation (5.12). □ 

Let us remark that in the homogeneous limit £j = 0, j = 1,...,N, and for two identical 
Bethe eigenstates, the evident mean value of can be easily derived from these representations. 
Indeed (5.12) and (3.19) yield 

n n 

<0| EI C{\ 3 ) n B(X k ) |0> 
q z M _ j=i k=i det(cD / ({A i })-2P / ({A J }) 

(0| n C{\j) n £( A <0 1°) detcp^A.ljj 
i=i fc=i 

where $'({Aj}) is the Gaudin matrix given by (3.20) and P'({Xj}) a matrix of rank one: 



Due to the following property of the Gaudin matrix, 

d 



6=1 

it can be rewritten as 

s z {n) = det n (J - 2C7), U ab = —, 

which leads to the evident value of s z (n) (U being a rank one matrix): 

, , s 2n 
^(n) = 1 - ^. 

5.3 The correlation function of spins at adjacent sites 

The next quantity to compute is the simplest two-point function, the correlator of two spins at 
adjacent sites, 

n n 

F- + (m,m + l,{^},{A fc })^(0| \{ CQij) a~ a+ +1 ]jB(X k ) |0). (5.18) 

j=l k=l 
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As usually {Afe} and {/ij} are supposed to be solutions of Bethe equations. From theorem 4.1 



this function can be written only in terms of the operators A, B, C and D. As previously, Bethe 
states being eigenstates for the propagator and the propagator from the first to the last site 
being equal to identity, we obtain a simple expression for this correlation function, 

F- + (m,m + 1, {fj,j}, {A fe }) = m _i({^}) <A m +i({Afe}) 

n n 

(0| l[C( H ) B(( m )C(Cm+i) l[B(X k )\0}. (5.19) 
j=i k=i 

It can then be reduced to a sum of scalar products by means of the commutation relations 
between the matrix elements of the monodromy matrix, namely, 

n n n n 

C(U+i)ilB(\ k )\0) = ^M a llB(\ k )\0) + ^M ab B(U+i) [J B ( X k)\Q)> 

k=l a=l k=l a ^b k=l (5.20) 

with the coefficients M a and M ab given by 

M - = n ^7 ~A X ^ II b-H^U+Jb-HK, A fe ), 

M ab = -— -d{\ a )b-\\ a ,\ b ) fl 6- 1 (A fc ,A i )6- 1 (A a ,A fc ). ( 5 - 21 ) 

V(Aa - gm+lJPKM - Sm+lJ A£ b 

Hence 

F-+(m,m + l,{nj},{X k }) = ^ m -i{{^j}) <P m l +1 {{\k}) x 

^ M a £ n (-{jUf}, {£ m , Ai, . . . , A a , . . . , A n }) + 

a=l 

+ MabSndfAj}, {Cm, 6n+l, Al , . . . , A a , Aft, . . . , A n }) j , (5.22) 
a+b ' 

where the hat means that the corresponding parameter is not present in the set. Since {[J>j} is 
a solution of Bethe equations, the scalar products S n ({fJ,j}, {Afe}) can be represented as deter- 
minants of size n according to theorem 3.1. Therefore, two-point functions at adjacent sites are 
obtained as a sum of determinants of size n. 



6 Conclusion 

In this article we have computed explicit determinant representations for form factors of the 
XXZ Heisenberg spin-^ inhomogeneous finite chain. These determinants are simply given in 
terms of usual functions of the parameters of the model. Moreover, an adjacent sites correlator 
has also been determined using similar techniques. The knowledge of the factorizing F-matrices 
from jl?]] was an essential ingredient to achieve this goal. It also sheds some new light on the 
algebraic structure underlying the Bethe Ansatz approach to correlation functions. In particular, 
multi-point correlators have been expressed in terms of expectation values (on the ferromagnetic 
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reference state) of operator entries of the quantum monodromy matrix: this was the result of 
explicitly solving the quantum inverse problem for local spin operators at any site of the chain. 

Let us stress also that this method allowed us to give a very direct proof of the scalar product 
formula between a Bethe eigenstate and an arbitrary state generated by the successive actions 
of the operators B. This formula is beautiful but very mysterious, since it involves the Jacobian 
of the eigenvalues of the transfer matrix with respect to the parameters of the Bethe states. 
Although the proof is very transparent, we do not know a satisfactory a priori explanation of it, 
and one feels that there should be a more direct understanding of this formula. 

What remains to be done, at least for the form factors we have computed, is to describe 
their thermodynamic limit and to compare the obtained results to those following the approach 
H . This will be done in a forthcoming publication, for the spontaneous magnetization of the 
XXZ chain. 

Another interesting question concerns the higher spin Heisenberg models. To deal with these 
more general cases, one would like to have at disposal the analogue of the factorizing i^-matrices, 
but here for the higher (fused) i?-matrices. This question is now under study. 
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Appendix A 

We give in this appendix the determinant representation for the partition function of the six- 
vertex model with domain wall boundary conditions, initially obtained in [40], and propose a 
direct proof for it, using the -F-basis. 

The partition function of the six-vertex model with domain wall boundary conditions corre- 
sponds to a special case of scalar product for the XXX or XXZ 1 /2 spin chain: 

JV N ( 



Zn ({A a }, {£j}) = | JJ ij j|-Bi..jv(Ai;£i, . . . ,£n) £i...iv(A2;£i) • • • 

•••/?;... v(Av: £,.... ,6v)}{n ti }• (A.I) 

It should be mentioned that the same partition function can be represented similarly as a matrix 
element of the products of operators C(A): 

Zn ({AqJ, {£?}) = | JJ Ij | |ci..jv(£i; Ai, . . . , Xn) Ci...n(&, Ai, . . . , Ajv) • • • 

...Ci.. J v(6v;A l! ... ,AAr)}{n^' }■ ( A - 2 ) 

Recursion relations for this function were obtained in fS8|| . It was shown that they define 
completely the function. The solution to these relations was given in [[5J as a determinant both 
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in the XXX and XXZ cases. For the XXX case one has the following result: 



nf =1 uLi (a. - j 

fc (£* ~ IL> ( 

where AA({A a }, {£•,}) is the N x N matrix given by 



z N ({a q }, fe}) = . m n a ^ detA/ "({ A «}' te» > ( A - 3 ) 



AT aj = ^ 77 — . ( A.4) 

In the XXZ case, the expression is similar: 

TT^-i TT^-i sinh(A a — £,•) 

^({AJ,{^}) = n 7 rr - ^7T r^detAP^ttAj,^}) , 

ni>fcSmh(C fc -Cj)n a >/3smh(A a -A^) 



where A/" xxz ({Aq,}, {£j}) is the N x N matrix given by 

sinh?] 

sinh(A Q - Cj + rj) sinh(A a - 



Kf Z = — — 7 =ct^ — ■ (A.6) 



We give now an explicit derivation of these representations, based on direct calculations in 
the .F-basis. Indeed, using the expression of operator B (or C), one obtains a new recursion 
formula for the partition function, which corresponds to a development of the determinant 
in (A. 3) or (A. 5). Explicit calculations are performed here in the XXX case, but they are quite 
similar in the XXZ case. 

More precisely, as the state { fljLi tj } (respectively { Yij=i lj }) * s invariant under the 
left-action of i*i.,jv (respectively the right-action of F^ 1 ^, the formula (A.l) can be directly 
written in the .F-basis: 

Zn = { II ^ ||-^l-iv(Ai;^i, • • • ,Cat) • • • #i..Jv(Ajv;£i, . . . ,6v) j j Y[ h j ' 

and, using the expression ( 2.29| ), we make 5i..jv(Ajv;£i, . . . , £jv) act on the state {njLi Tj'}j i n 
order to obtain a recursion relation for 

z N =j2 c (a^&) (n & (a^v> go) { n is } {bi...n(*i) . ■ ■ bi.m*n-i)} { ( n ti ) uo } • 

We thus extract, for each term i of the sum, the action (which is now diagonal for (o"j~) 2 = 0) 
on space i of the other operators -Bi...jv(A q ), 1 < a < N — 1, which leads to an extra numerical 
factor: 



t=i S^i y 

I II ^ ||-Ei-i-ii+i...;v(Ai;£i, • • • • • • >£iv) • • • 

• • • £i...i_ij + i..jv(Ajv-i;£i, • • • • • • j6v) | j Y[ 1" 
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Indeed, as the number of operators B was formerly equal to the number of sites on the chain, 
and as (&J) 2 = 0, the product IljLi a J appears in all the non-zero terms in the development of 
the product of B. 

We have obtained the following recursion formula for Zn'- 

N , N v 

i=i S =a J (A.7) 

This corresponds actually to the last line development of the determinant in the formula 

7 !t\ 1 ft l\ nf=llla=l(^at - £j) y 
Zn ({A q }, {Cj}) = i=f 77 7VFF 7\ TT det -™ 



where AT is the matrix obtained from J\f by adding to the last line Ln the linear combination 

N . j. N-l 



of the other lines Y1b=i fpLp, with coefficients 



Ajv — A Q 



fe =i Ajv ~ ^ k + ^ «=i X P ~ Xa 
In fact, this development leads to the following recursion formula for Zn' 

z ^S^(S^)(n^)(n^|){n(v- &+ ,) 

k^i k^i k^i 

~ Yl fl( X P-^k+r]) \ N _ f f Yl \ N _\ a ) \ Z N -i{{K} a ^N,{ij}j^i) » 

0=1 fc=l ^ ?i \ a=l A P Aa ) J 

and one can easily see that 



N /N-l j. . \ ( N 



nK.-^+")=(nfe^]m^-& + ") 



N-l N , A /JV-1 



equality between two polynomials of degree iV — 1 in £j, which can be proved at the N points 
& = A a , 1 < a < iV. 

Thus, as A/" et jv have the same determinant, this finishes our proof of the formula ( |A.3| ) for 
the partition function. 
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Appendix B 

For the calculation of scalar products and correlation functions in the -F-basis, we need deter- 
minant representations for the following functions 



G^ ) ({X k },i 1 ,...,i n ) = (i 1 ,...,i n \ JjB(A fc )|0), (B.l) 

k=l 

n 

GP({ fM },i 1 ,...,i n ) = (0\ l[C(jn) \h,...,i n )- (B.2) 



Such representations can be easily obtained from the one of the partition function Z n , by cal- 
culating explicitly the action of the operators Si...jv(A a ) or Ci..jv(^ a ) in the sites which do not 
belong to . . . ,i n }. Indeed, in the expression for the operator B(X) in the F-basis 



N 

~[ m V b («»«) J [k] 

only the terms corresponding to a cr~ , where k belongs to {1, . . . ,n}, give a non-zero contribution 

to the function Gg ({A^}, i±, . . . ,i n ). It means that the operators Bi j^(X a ), 1 < a < n, 
act as diagonal ones on all the spaces except ii, ■ ■ ■ ,i n - Thus, we can extract the action of 
n«=i ^l...Ar(A a ) on all the sites but i\, . . . , i n as a global factor: 

/ n N 



d§ ({A fc }, h, . . . , i n ) = J] If b(\ a ,£ k ) Z n ({A Q }, }) 



ct=l fc=l 

n n 

n n N n n^^-^+i) 

IIII *(*«>&)) n "7^"% S n n rrdet^({A a },fe,}) , 

t=i J n y - Cij n ^ (Aa - a^) (B 4) 



where A/"({A Q }, is the n x n matrix defined by (A. 4) or (A. 6). 

Using t 
representa! 
is given by 



Using the representation of the operator C(/i) in the F-basis one can analogously obtain a 

(o) 
c 



representation for the function Gq ({/i/}, i n ) (B.2). Both in the XXX and XXZ cases it 



/ n N \ ( n \ 

Gg ) ({ W },i 1 ,...,i n )= m n 6"%,^) n )Gf(h,...,i n ,{»i})- 

S=l fe=i ' S\fc=i ' (B.5) 

Appendix C 

In this Appendix we calculate recursively the intermediate functions G ^ (|3.9| ). The detailed 
proof is written here only for the XXX case, but is quite similar in the XXZ case. 

To compute the function G W starting from the expression ( 3.12j ) for G ^ , one has to perform 



the summation in (3.10) . Taking into account that only the first column of the matrix J\f depends 
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on the parameter ^ , one can express G W as follows 

G ( - 1) ({X a },fi 1 ,i 2 , ...,i n ) = — — — Z~ lk ~ 2 ^ \ 7T — det7V (1) ({A a },^i,i 2 , • • • ,i n ), 

11 (« fc -«J 11 (A/J-Aa) / C1 s 
n>j>k>2 l<a<f3<n V " ' 

where the matrix TV^^jAa}, //i, Z2, . . . ,i n ) is defined as 

^ab=—^Mab for 6>2, (C.2) 



n n 

fc=2 

It should be mentioned that the summation in (C.3) can be taken over all the possible values of 
i\ as the contributions of the terms i\ = ij, j > 1 are equal to zero. It is possible to calculate 
explicitly the sum in (C.3) using its analytical properties and the fact that {X a } is a solution of 
Bethe equations: 

n (Am - & +v) 

E V V ™¥=a T-r h -\ (t . f \ tigx 



1 ft(Am-eo 

1 m=l '/ 



where the function i? a & has the form 

Hab = T~Z — ( r (^ 6 ) Y\( Xm ~ Vb + V) ~ n (A m - //ft - */) ) , (C- 5 ) 

Aa A*6 ^ , , ' 

m^a m^=a 

with 

Indeed, the left hand side of (C.4) is a rational function of \x\ with simple poles at the points 
Mi = £jj i = Ij ■ ■ ■ > N, and its limit is zero when \i\ — > oo. The right hand side is also a rational 
function of /xi, which has only simple poles and becomes zero when /xi — > oo. The residues of 
the r.h.s. at the points = £j are the same as in the l.h.s. Thus one should only prove that 
the r.h.s. has no other poles, namely when \i\ = X a and /xi = £j fc — rj, k = 2, . . . , n. One can 
easily see that the residues of the r.h.s. at the points \i\ = £i k — V are equal to zero. As {A Q } is 
a solution of Bethe equations, the residue at the point [i\ = \ a is also equal to zero. Therefore 
the l.h.s and r.h.s of (C.4) are rational functions having the same behavior when fi\ — > oo, the 
same simple poles and the same residues in these poles, thus they are equal. 
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So, the matrix elements of the first column of the matrix AA 1 ) have the form 

6=2 



where a b are coefficients which do not depend on a. Only the first term in this sum gives a 
nonzero contribution to the determinant of J\f^ , which leads to the following representation for 



n n 



n n (Aa-^+r?) 

G {1) ({X a },ni,i2, ... , i„) = — — " fc ~ 2 ^ \ 7T — det£ (1) ({A Q },^i,i 2 ,. . . ,i„), 

11 U>. i'J 11 (A/3-A a J / Cfi x 

n>j>k>2 l<a</3<n K^-^J 

with the matrix C/^ 1 ) defined as 



Mi " & (C.7) 



Repeating this procedure, one obtains a general expression for all the functions G ( m ' : 

n n 

n n (A Q -e, fc + ^) 



<u> U A fe/! Ml) ■ ■ ■ i Mmj^m+l) = - 



n o,) n (A/s-Aa) 

n>j>A:>m+l l<Q</3<n 

dett/ (m) ({A Q },/xi, . . . ,/i m ,i m+ i, . . . ,i n ), (Ci 



n (pk - h) 

m>j>k>l 



with the matrix £( m) ({A Q }, /ii, . . . ,/x m ,«m+i, . . . ,i n ) given by 

slr^^n H for 6 >™> ( c - 9 ) 

G^=H ab for 6<m. (CIO) 



We prove (C.8)-(C.10) by induction. For m = 1 it coincides with ( p.6| ). Let this representa- 
tion be valid for G (m_1) . Combining it with (3.10) and (gig ) we obtain the following expression 
for G^: 



n n 



n n {K-^ k +v) 

^m+l j • • • j In) 



tm)(f \ \ ■ ■ \ _ q=l fc=m+l 

ll^fc/j Ml) — j mi *to+1j — j In) — -p-r j. \ -p-r \ \ ' N 

11 \$ik-&j) 11 (A/3-A Q J 

n>j>fc>m+l l<a<(3<n 

det7V (m) ({A Q },^i, . . . ,fjL m ,i m+ x, . . .,i n ), (C.ll) 



n (m* - Mi) 

m— l>j>fc>l 
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where the matrix M^ m \{X a }, fix, . . . , fi m , i m +x, . . . ,i n ) is defined as 



for b > m, 



Mr ] =H ab for b<m, 



3 

'ah = n oi> 

n n 

*ffl= II (Mi-^+»y)E ; h x h x 

fc=m+l i^l^" 1 ^ mAa ^ m 



n - &» + v) 

x — ~i II b-HUUi)- (C12) 

n - + »7) n - 

fc=m+l J=l 

The sum in (C.12) can be computed the same way as in (C.4). One can prove using similar 
arguments that 

„ m—1 n 

\f(m) _ f*am , \ ~* R (m) u \ ^ (m) . r 

n(^ -Mm) n (v m -Zi k +v) b=i b = m+i ( c - 13 ) 

3=1 fc=m+l 

where a^ 1 and /Sj' are coefficients which do not depend on a. As only the first term gives a 
nonzero contribution to the determinant we obtain the representation (C.8). 
Finally the scalar product is given by 

Sn({M,{A a }) = GW({A a },/, 1 ,..., Mn )- detH({X a },{ H }) 



n (m* - Mi) n (A/? - Aa) ' (c u) 



j>k a<f3 

with matrix elements of H({X a }, {/x?}) defined by ( |C . 5|) . 
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